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We  report  thermal  rate  coefficients  for  12  reactions  of  rare  gas  cations  (Ne+,  Ar+,  Kr+,  Xe+)  with 
halide  anions  (Cl”,  Br”,  I”),  comprising  both  mutual  neutralization  (MN)  and  transfer  ionization. 
No  rate  coefficients  have  been  previously  reported  for  these  reactions;  however,  the  development  of 
the  Variable  Electron  and  Neutral  Density  Attachment  Mass  Spectrometry  technique  makes  it  pos¬ 
sible  to  measure  the  difference  of  the  rate  coefficients  for  pairs  of  parallel  reactions  in  a  Flowing 
Afterglow-Langmuir  Probe  apparatus.  Measurements  of  18  such  combinations  of  competing  reac¬ 
tion  pairs  yield  an  over-determined  data  set  from  which  a  consistent  set  of  rate  coefficients  of  the 
12  MN  reactions  can  be  deduced.  Unlike  rate  coefficients  of  MN  reactions  involving  at  least  one 
polyatomic  ion,  which  vary  by  at  most  a  factor  of  ~3,  those  of  the  atom-atom  reactions  vary  by  at 
least  a  factor  60  depending  on  the  species.  It  is  found  that  the  rate  coefficients  involving  light  rare- 
gas  ions  are  larger  than  those  for  the  heavier  rare-gas  ions,  but  the  opposite  trend  is  observed  in  the 
progression  from  Cl”  to  I”.  The  largest  rate  coefficient  is  6.5  x  10”^  cm^  s”^  for  Ne+  with  I”.  Rate 
coefficients  for  Ar+,  &■*■,  and  Xe+  reacting  with  Bra”  are  also  reported.  ©  2014  AIP  Publishing 
LLC.  [http;//dx.doi.org/10. 1063/1. 4862151] 


INTRODUCTION 

The  mutual  destruction  of  positive  ions  and  negative  ions 
in  ion-ion  plasmas  can  occur  by  several  mechanisms.  Ter- 
molecular  ion-ion  recombination,  in  which  energy  is  trans¬ 
ferred  to  ambient  gas  atoms  M, 

A+ +  P- +  M  ^  A  + B  +  M,  (la) 

-^AB  +  M  (lb) 

has  been  treated  theoretically  in  great  detail  by  classical  me¬ 
chanics  (see,  e.g.,  Bates  or  Flannery).'’^  In  that  case,  the 
chemical  nature  of  the  recombining  ions  enters  only  weakly 
through  their  masses  and  ion-atom  interaction  potentials.  In 
contrast,  binary  mutual  neutralization  (MN), 

A+  +  B-  ^  A  +  B  (2) 

depends  sensitively  on  the  internal  and  electronic  states  of 
the  interacting  species.  When  both  ions  are  atomic,  MN  must 
occur  at  intersections  of  the  initial  ionic  potential  curve  and 
those  describing  the  interaction  between  the  atoms  in  their  fi¬ 
nal  states,  as  in  “accidentally  resonant”  charge  transfer  of  ions 
and  atoms.  Unlike  reactions  involving  a  polyatomic  species, 
which  will  have  a  broad  manifold  of  possible  product  elec¬ 
tronic  and  vibrational  states,  an  atom-atom  system  will  have 
a  sparse  density  of  possible  product  electronic  states.  The 
rate  coefficients  for  MN  then  depend  critically  on  the  chemi¬ 
cal  species  involved  and  are  not  easily  predictable  by  simple 
arguments. 
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With  the  exception  of  the  H+  -|-  H”  system,  for  which 
cross-sections  have  been  determined  by  experiment^  and 
calculation,^  and  a  handful  of  other  systems,  experimental 
data  on  binary  atomic  MN  reactions  are  rather  sparse.  Rate 
coefficients  of  atomic  rare  gas  cations  with  halide  anions 
(Rg+  +  X”)  reactions  have  been  measured  in  high-pressure 
afterglows^  but  the  extrapolation  to  low  pressures  is  highly 
uncertain.  There  are  quite  a  few  theoretical  estimates  based 
on  the  Fandau-Zener  model,®’ ^  and  upper  limits  on  four  bi¬ 
nary  Rg+  +  X”  MN  rate  coefficients  have  been  determined 
experimentally  in  a  flowing  afterglow.*  Accurate  theoretical 
treatments  are  inherently  difficult  because  curve  crossings  in¬ 
volving  excited  state  species  at  large  internuclear  distance  can 
contribute.  Approximate  scaling  formulas^’ are  useful  to  ob¬ 
tain  rough  estimates,  but  they  rely  on  simplifying  assump¬ 
tions  that  may  not  be  justified  in  all  cases.  The  Double  Elec¬ 
trostatic  Storage  Ring  Ion  Experiment  will  soon  be  brought 
online"  and  is  expected  to  yield  important  new  data  on  mu¬ 
tual  neutralization;  much  like  the  advent  of  magnetic  storage 
rings  in  the  1990s  greatly  advanced  the  study  of  dissociative 
recombination. 

Recently,  we  have  used  the  variable  electron  and  neu¬ 
tral  density  attachment  mass  spectrometry  (VENDAMS) 
technique"  to  study  a  wide  variety  of  MN  processes.®  These 
include  the  first  complete  product  distributions"  and  evidence 
that  electrons  catalyze  the  process."’"  The  measurements 
rely  on  the  difference  in  rate  coefficients  of,  typically,  Ar+ 
with  the  anion  of  interest  to  that  of  Ar"*"  with  an  atomic  halide 
anion,  the  latter  of  which  have  been  assumed  to  proceed  at 
a  negligible  rate.*’®  In  this  paper,  using  a  slightly  modified 
VENDAMS  method,  we  describe  the  results  of  a  series  of 
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measurements  on  binary  mutual  neutralization  reactions  of 
atomic  rare  gas  cations  with  atomic  halide  anions. 

Direct  measurement  of  either  cation  or  anion  absolute 
densities  under  typical  conditions  in  a  flowing  afterglow- 
Langmuir  probe  (FALP)  apparatus  using  the  standard 
“orbital-motion”  theory  is  problematic  due  to  collisions 
between  ions  and  neutral  species  trapping  ions  near  the 
probe  and  thereby  artificially  increasing  the  probe’s  target 
area.'®  Additionally,  measurement  of  any  MN  rate  coeffi¬ 
cients  smaller  than  ~10“*  cm^  s“*  is  difficult  in  traditional 
FALP  measurements  due  to  the  competing  rate  of  ambipo- 
lar  diffusive  loss.  However,  measurement  of  the  relative  con¬ 
centration  of  either  anions  or  cations  in  FALP  mass  spectra 
is  routine  and  robust.  By  creating  conditions  with  a  single 
dominant  cation  and  two  or  more  anionic  species,  the  change 
in  the  relative  abundance  of  the  anions  as  a  function  of  the 
initial  plasma  density,  with  which  the  initial  cation  density 
scales,  indicates  the  difference  in  the  rate  coefficients  of  the 
competing  neutralization  reactions.  Importantly,  only  electron 
densities  need  to  be  measured  absolutely,  with  all  ion  density 
measurements  being  relative.  Similar  experiments  can  be  con¬ 
ducted  by  creating  conditions  with  a  single  dominant  anion 
and  multiple  cations.  By  making  difference  measurements  of 
many  combinations  of  both  anions  (Cl~,  Br“,  I~)  and  cations 
(Ne+,  Ar+,  Kr+,  Xe+),  a  ladder  of  reactivity  for  the  series 
of  reactions  Rg+  -|-  X“  is  established.  Because  the  same  re¬ 
actions  are  probed  via  monitoring  either  competing  cations 
or  competing  anions,  redundancy  results  in  the  ladder  being 
over-determined,  with  internal  consistency  giving  increased 
confidence  in  the  resulting  values.  Third-body  contributions 
due  to  neutral  atoms  are  negligible  at  the  experimental  gas 
densities  (on  the  order  of  3  x  10*®  cm^  s“'  of  He),  as  shown 
in  a  recent  systematic  study  of  a  possible  enhancement  of  MN 
rates. *^  Similarly,  both  the  densities  of  ions  (on  the  order  of 
10® — 10*®  cm^  s“*)  and  electrons  are  too  small  to  affect  the 
rate  coefficients. 

In  some  of  the  reactions,  the  halogen  product  is  formed 
in  excited  states  with  sufficient  energy  to  autoionize  to  form 
X+  -Fe", 

Rg+  +  X-  ^  Rg  +  X*,  (3a) 

^Rg  +  X+  +  e-.  (3b) 

This  process  (3b)  has  not  been  reported  previously  for 
cation-anion  neutralization,  but  is  similar  to  transfer  ioniza¬ 
tion  (Tl)  observed  in  other  classes  of  charge  transfer.  Results 
for  TI  of  Ar+  to  Br~  and  H  are  reported  in  Ref.  18.  The 
rate  coefficients  reported  here  are  for  the  sum  of  MN  and  TI, 
i.e.,  reactions  (2),  (3a),  and  (3b)  although  Tl  is  endoergic  and 
cannot  occur  in  about  half  of  the  reported  systems. 

EXPERIMENTAL 

The  measurements  were  made  using  the  Air  Force 
Research  Laboratory’s  FALP  apparatus  using  the  VEN- 
DAMS  technique,  both  of  which  have  been  described 
in  detail  previously.*^  *®  At  the  upstream  end  of  the  1 
m  long  FALP,  a  helium  (99.999%,  Matheson)  flow  of 
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1 0-20  std.  L  min“  *  is  ionized  by  a  microwave  discharge  yield¬ 
ing  He+,  He2''',  He*,  and  e“.  An  Ar  (99.999%,  Matheson) 
flow  of  ~0.3  sfd.  L  min“*  is  added  ~10  cm  downsfream  con¬ 
verting  He2'''  and  He*  to  Ar+,  yielding  a  plasma  in  which  over 
95%  of  the  cations  are  Ar+  with  the  remainder  being  He+  and 
small  levels  of  water  and  air  ion  impurities.  By  replacing  the 
Ar  with  Rr  or  Xe,  a  similarly  pure  Kr+  or  Xe+  plasma  can 
be  made.  The  reaction  of  He*  with  Ne  to  yield  Ne+  is  en¬ 
dothermic,  precluding  production  of  a  Ne+  dominated  plasma 
in  the  same  manner.  Alternatively,  two  rare  gases  instead  of 
one  may  be  added  downstream,  each  through  a  separate  mass 
flow  meter  but  mixing  prior  to  entering  the  flow  tube  through 
a  single  inlet,  to  yield  a  mixture  of  Rg+  cations.  Experiments 
were  conducted  in  Ar+/Ne+,  Ar+/Xe+,  and  Ar+/Kr+  plas¬ 
mas  with  typically  >80%  Ar+,  15%  Ne+/Xe+/Kr+,  and  <5% 
other  cations.  Typical  EALP  studies  are  conducted  at  ~  1  Torr 
at  300  K;  this  study  was  conducted  at  higher  pressure  (1.5-2 
Torr)  in  order  to  convert  a  larger  fraction  of  He+  into  He2’*", 
which  then  converts  rapidly  to  Rg®^,  thereby  reducing  the 
fraction  of  He+  in  the  plasma. 

Either  one  or  two  neutral  attaching  gases  (CCI4  (Sigma 
Aldrich)  to  yield  CD,  CE2Br2  (Sigma  Aldrich)  to  yield  Br^, 
or  C2E3I  (SynQuest  Labs)  to  yield  I“)20'2i  may  be  added  to 
the  afterglow  through  a  glass  inlet  comprising  6  hngers  evenly 
spaced  on  a  1  cm  radius  located  44  cm  before  the  end  of 
the  flow  lube.  These  reagents  are  prepared  as  dilute  (0.5%- 
1.0%)  mixtures  in  helium  and  added  using  mass  flow  me¬ 
ters.  The  liquid  reagents,  CCI4  and  CF2Br2,  are  purihed  via 
freeze-pump-thaw  cycles  prior  to  mixture  preparation. 

The  absolute  electron  density  along  the  center  axis  of 
the  flow  tube  may  be  measured  using  a  cylindrical  Lang¬ 
muir  probe.  The  probe  is  movable  from  15  cm  before  to 
35  cm  after  the  neutral  inlet;  however,  with  the  exception  of 
diffusion  measurements,  all  electron  density  measurements 
here  are  made  1  cm  upstream  of  the  neutral  reactant  inlet, 
and  reflect  the  initial  electron  density  [e]o  upon  introduc¬ 
tion  of  reactant  species.  [e]o  may  be  varied  over  several  or¬ 
ders  of  magnitude  (in  the  present  case  from  ~10*  cm“®  to 
~2  X  10*®  cm^®)  primarily  by  translating  the  microwave  dis¬ 
charge  cavity  farther  from  or  closer  to  the  reactant  inlet.  Im¬ 
portantly,  and  in  contrast  to  prior  EALP  studies  of  MN,*  the 
Langmuir  probe  is  used  here  only  to  measure  electron,  and 
not  ion,  densities.  Prior  to  the  introduction  of  any  electron  at¬ 
taching  species,  the  total  cation  concentration  is  assumed  to 
be  equal  to  the  electron  concentration  (i.e.,  the  plasma  is  neu¬ 
tral),  and  [e]o  is  taken  as  also  a  measurement  of  the  plasma 
density.  The  center  axis  of  the  flow  al  the  downstream  end  of 
the  flow  tube  is  sampled  through  a  0.3  mm  diameter  aperture 
in  a  truncated  nosecone  leading  to  a  lens  stack  and  quadrupole 
mass  spectrometer.  Ion  velocity  and  ambipolar  diffusion  are 
measured  using  standard  methods.  Diffusion  after  conditions 
evolve  into  an  ion/ion  plasma  is  assumed  to  occur  with  a  rate 
coefficient  one  half  that  measured  in  the  cation/e“  plasma.22 
Mass  discrimination  is  not  particularly  important  to  the  cur¬ 
rent  experiments,  but  is  measured  where  needed  by  methods 
previously  described. 2* 

Measurements  were  made  for  each  pair  of  Cl“ ,  Br“ ,  and 
1“  in  each  of  the  plasmas  dominated  by  either  Ar+,  Kr+,  or 
Xe+.  Experiments  were  attempted  using  0.5%  F2  in  He  as 
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a  source  of  F“;  however,  significant  chemistry  appeared  to 
occur  between  F2  and  the  other  neutral  precursors  in  the  in¬ 
let  lines,  resulting  in  unknown  neutral  concentrations  enter¬ 
ing  the  flow  tube  and  precluding  an  accurate  analysis  of  the 
data.  Measurements  were  also  made  in  Ar+/Xe+,  Ar+/Kr+, 
and  Ar+/Ne+  with  each  of  the  anions  Cl“,  Br” ,  and  1“ .  These 
experiments  were  not  attempted  with  F“ ,  as  F2  attaches  elec¬ 
trons  inefficiently  preventing  the  quick  depletion  of  electron 
density  with  a  low  concentration  of  added  neutral;  this  is  an  is¬ 
sue  as  the  sampling  efficiency  and  mass  discrimination  of  the 
detection  system  changes  significantly  between  a  plasma  with 
considerable  concentrations  of  electrons  and  one  without.  The 
former  would  occur  at  high  [e]o  and  the  latter  at  low  [e]o- 
Nascent  ratios  of  Ar+;Rg+  as  a  function  of  [e]o  were  mea¬ 
sured  without  any  attaching  gas  added.  For  both  the  Ar+/Xe+ 
and  Ar+/Kr+  plasmas,  this  ratio  was  independent  of  [e]o.  The 
ratio  in  the  Ar+/Ne+  showed  structure  with  [e]o  as  discussed 
below. 

DATA  AND  ANALYSIS 

The  primary  data  in  VENDAMS  are  relative  ion  concen¬ 
trations  measured  after  a  fixed  reaction  time,  typically  7.1  ms, 
as  a  function  of  [e]o,  as  measured  using  the  Langmuir  probe 
just  before  the  neutral  inlet.  The  present  measurements  were 
performed  one  of  two  ways,  either  monitoring  the  concen¬ 
trations  of  two  cationic  species  in  the  presence  of  a  single 
dominant  anion,  or  the  inverse,  monitoring  the  concentrations 
of  two  anionic  species  in  the  presence  of  a  single  dominant 
cation. 

For  the  former  technique,  two  noble  gases  (e.g.,  Ar  and 
Kr)  are  added  to  an  upstream  inlet  to  create  an  Ar+/Kr+/e“ 
plasma  as  described  above.  Downstream,  a  single  gas  that 
rapidly  attaches  electrons  is  added  in  known  concentration 
through  the  neutral  inlet  using  a  mass  flow  meter  to  convert 
all  electrons  into  halide  anions  within  a  few  cm;  the  complete 
electron  depletion  is  confirmed  using  the  Langmuir  probe. 
CCI4,  CF2Br2,  and  C2F3I  were  used  to  make  CD,  Br~,  and  D, 
respectively. The  neutral  concentration  is  kept  sufficiently 
low  (typically  ~3  x  10^*'  cm“^)  that  reaction  with  Rg+  occurs 
to  only  a  small  extent,  and  the  cationic  species  in  the  plasma 
remain  predominantly  Rg’*’.  Representative  data  for  Ar+/Kr+ 
reacting  with  1“  are  shown  in  Fig.  1.  As  the  plasma  density 
increases,  the  absolute  concentration  of  1“  formed  via  attach¬ 
ment  increases  and  the  abundance  of  Ar+  decreases  relative  to 
that  of  Kr+ ,  reflecting  a  larger  rate  of  reaction  between  Ar+ 
and  D  than  between  Kr+  and  D .  In  the  particular  case  of  the 
Ar+/Ne+  plasma,  the  nascent  fraction  of  Ne+  in  the  plasma 
varied  with  [e]o  (Fig.  2).  This  background  variation  was  in¬ 
corporated  into  the  modeling  described  below,  increasing  the 
uncertainty  in  the  derived  values. 

The  analogous  anion  measurements  involve  instead  a 
plasma  dominated  by  a  single  cation,  either  Ar+,  Xe+,  or 
Kr+  as  described  above.  Two  attaching  gases  are  added  in 
known  concentrations  through  separate  mass  flow  meters,  but 
through  the  same  neutral  inlet,  yielding  two  atomic  halide  an¬ 
ions.  The  relative  concentrations  of  the  anions  are  monitored 
as  a  function  of  [e]o,  yielding  information  on  the  reactions 
of  those  anions  with  the  dominant  cation.  Representative  data 
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FIG.  1 .  Relative  cation  abundances  as  a  function  of  initial  electron  density 
7.1-ms  after  addition  of  6.5  x  10^^  cm~^  C2F3I  to  an  Ar'^/Kr'^/e  plasma  at 
300  K.  Solid  lines  are  best  fit  models;  dashed  lines  are  representative  fits  at 
the  reported  uncertainty  limits. 

from  the  addition  of  CCI4  and  CF2Br2  to  an  Ar+/e“  plasma 
is  shown  in  Fig.  3.  The  Br2“  is  due  to  a  7%  branching  in  the 
electron  attachment  to  CF2Br2;^°  however,  the  presence  of  the 
minor  competing  ion  does  not  add  significant  difficulty  to  the 
analysis,  and  in  fact  adds  ancillary  information  on  the  reac¬ 
tion  of  Br2~  with  the  Rg+  species.  Interpretation  of  the  data 
here  is  similar  to  the  cation  case,  indicating  that  Br“  and  CD 
react  with  Ar+  at  similar  rates,  with  Br”  being  slightly  faster, 
and  Br2“  reacting  significantly  more  rapidly  with  Ar+. 

While  there  is  no  simple  analytic  form  to  extract  rate 
coefficients  from  the  data,  information  on  the  rate  coefficients 


FIG.  2.  Relative  cation  abundances  as  a  function  of  initial  electron  density 
7.1  ms  after  addition  of  3.5  x  10^®  cm~^  C2F3I  to  an  Ar+/Ne'^/e  plasma 
at  300  K  (solid  points)  and  under  the  same  conditions  without  addition  of 
C2F3I  (open  points).  Solid  lines  are  best  fit  models  accounting  for  the  nascent 
distribution  of  Ar'^iNe'*';  dashed  lines  are  representative  fits  at  the  reported 
uncertainty  limits. 
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FIG.  3.  Relative  anion  abundances  as  a  function  of  initial  electron  density 
7.1  ms  after  addition  of  9.3  x  10®  cm~^  CCI4  and  2.2  x  10*®  cm“^  CF2Br2 
to  an  Ar'*'/e  plasma  at  300  K.  Solid  lines  are  best  fit  models;  dashed  lines  are 
representative  fits  at  the  reported  uncertainty  limits. 

of  the  competing  reactions  may  be  derived  by  numerically 
solving  the  rate  equations  describing  the  chemistry  occurring 
in  the  flow  tube  from  the  known  initial  conditions  throughout 
the  known  reaction  time.  As  described  in  detail  in  Ref.  12, 
the  full  parameter  space  defined  by  uncertainties  in  rate  co¬ 
efficients,  neutral  concentrations,  and  diffusion  coefficients  is 
sampled  through  a  Monte  Carlo  method,  with  calculated  rel¬ 
ative  ion  abundances  compared  to  the  measured  data  through 
a  weighted  least  squares  fit.  While  the  complete  chemistry 
is  modeled,  the  data  are  almost  exclusively  sensitive  to  dif¬ 
ference  of  the  rate  coefficients  (comprising  both  MN  and  TI 
product  channels)  for  the  two  competing  reactions  and  the 
measured  diffusion  rate.  The  latter  affects  the  total  concen¬ 
tration  of  the  species  and  therefore  the  rate  of  disappearance. 
It  should  be  noted  that  all  ions  were  treated  as  diffusing  at  the 
same  rate.  It  is  unclear  if  this  assumption  is  justified;  however, 
that  the  results  of  these  measurements,  which  involve  differ¬ 
ent  pairs  of  ions  yielding  information  on  the  same  reaction 
rate  constants,  were  self-consistent  suggests  that  differential 
diffusion  rates  between  the  ions  are  small. 

The  data  from  each  individual  experiment  here  constrain 
the  difference  in  the  competing  rate  coefficients,  as  opposed 
to  their  absolute  values.  These  results  are  compiled  in 
Tables  I  and  II.  Because  there  is  overlap  in  the  reactions 
studied  between  experiments,  these  differences  form  a  ladder 
connecting  the  set  of  Rg+/X“  rate  coefficients.  Because  of  the 
redundancy  in  the  measurements  between  the  two  techniques 
(monitoring  either  anions  or  cations),  the  interconnected  set 
is  over-determined  (Fig.  4(b)).  That  the  over-determined  set 
is  self-consistent  supports  the  validity  of  the  measurements. 

If  any  one  of  the  rate  coefficients  in  the  ladder  is  known, 
then  all  the  values  can  be  put  on  an  absolute  scale.  No  two- 
body  rate  coefficients  for  Rg'''/X“  systems  have  been  reported 
in  the  literature;  however.  Church  and  Smith  report  upper 
limits  of  0. 2-1.0  x  10“*  cm^  s“'  for  reactions  of  F“  with 
Ar+,  Kr+,  and  Xe+,  and  an  upper  limit  of  5  x  10“®  cm^  s“' 
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TABLE  I.  Measured  differences  in  the  rate  coefficients  of  Ar'®  -|-  X  and 
Rg+  -bX-. 


X- 

Rg+ 

^Ar++x-  -  kRg++x-  (10  *  s  *) 

cr 

Ne+ 

-2.6  ±  1.1 

Kr+ 

U.DUlo  1 

Xe+ 

U.^ili02 

Br- 

Ne+ 

-2.8  ±  1.2 

Kr+ 

U.DIII02 

Xe+ 

U.DUlos 

I 

Ne+ 

-4.5  ±  1.4 

Kr+ 

1.2  ±0.2 

Xe+ 

2 

^•^=^0.4 

for  Xe+  -|-  Cl“ ,  the  latter  of  which  appears  in  our  data  set  pro¬ 
viding  a  constraint  for  the  present  results.^  Further,  our  ladder 
limits  the  rate  coefficients  for  several  of  the  reactions  (Kr+ 
-|-  Cl“,  Kr+  +  Br“,  and  all  three  Xe+  reactions)  to  be  the 
same  value  within  about  2  x  10“^  cm^  s“',  the  sensitivity  of 
the  experiment.  While  it  is  possible  that  all  of  these  reactions 
occur  with  similar  rate  coefficients  close  to  the  reported  upper 
limit  for  Xe+  -|-  CD,  it  appears  more  likely  that  all  the  reac¬ 
tions  simply  occur  with  a  rate  coefficient  below  the  detection 
limit  of  the  experiment.  If  we  assume  that  any  one  of  these 
five  reactions  occurs  with  a  rate  coefficient  below  ~10“^  cm^ 
s“',  then  this  becomes  an  effective  anchor  for  all  the  mea¬ 
sured  rate  coefficients  as  it  is  much  less  than  the  experimental 
uncertainty.  Therefore,  in  addition  to  the  18  analyses  of  the 
individual  data  sets  yielding  the  difference  measurements 
reported  in  Tables  I  and  II,  we  have  also  analyzed  15  data 
sets  (excluding  those  involving  Ne"*",  which  lack  redundancy) 
in  concert  with  the  constraint  that  Xe+  -|-  Cl“  proceeds  with 
k  <  10“®  cm^  s“^.  The  derived  rate  coefficients  are  shown  in 
Fig.  4(a)  and  listed  in  Table  III.  The  set  of  measurements  was 
self-consistent,  as  demonstrated  in  Fig.  4(b)  comparing  the 
individual  difference  measurements  to  the  derived  rate  coef¬ 
ficients;  an  inconsistent  data  set  plotted  in  this  manner  would 
show  some  of  the  termini  of  the  ball  and  chains  outside  of  the 
shaded  uncertainties  for  the  derived  absolute  rate  coefficients. 
It  is  worth  noting  that  the  Ne+  rate  coefficients  are  each 


TABLE  II.  Measured  differences  in  the  rate  coefficients  of  the  indicated 
Rg'^  +  X“  and  Rg”^  +  Y“  reactions. 


^Rg+AX- 

—  ^Rg++Y-(^^  * 

cm^  s 

X- 

Rg+ 

Br- 

r 

Br2" 

Ar+ 

cr 

^•■^='=0.3 

3.4  ±  1.0 

Br- 

^•^='=0.3 

3.2  ±  0.9 

I- 

1.7  ±0.6 

Kr+ 

cr 

-0.1  ±  0.4 

o.8±[>:« 

3.3  ±0.8 

Br- 

0.9±[>;« 

3.2  ±0.8 

I- 

2.3  ±  0.8 

Xe+ 

cr 

0  2±^'^ 

0.1  ±0.2 

4.9  ±  1.6 

Br- 

0-2±'i;? 

4.0  ±  1.3 

I- 

3.0  ±  1.6 
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FIG.  4.  (a)  Derived  rate  coefficients  for  reactions  of  the  indicated  Kg'**  +  at  300  K.  Triangles  indicate  upper  limits,  (b)  Representation  of  the  over-determined 
ladder  of  measured  differences  in  the  Rg''*  -|-  X~  rate  coefficients  (see  text).  The  derived  rate  coefficients  are  shown  as  in  (a),  and  are  overlayed  with  the  measured 
differences  (ball  and  chains)  between  the  rate  coefficients  of  the  indicated  reactions.  Color  of  the  balls  indicates  the  relevant  anion  (green  =  CD;  brown  = 
Br“;  black  =  D),  while  the  shaded  regions  at  each  terminus  indicate  the  cation  (red  =  Ne+;  purple  =  Ar*^;  blue  =  Kr+;  grey  =  Xe+).  The  extents  of  the 
chains  are  the  measured  differences  (typical  uncertainty  is  20%,  not  shown).  Vertical  extents  of  the  shaded  areas  are  at  the  uncertainty  limits  for  each  best-fit 
rate  coefficient;  horizontal  extents  of  the  shaded  areas  are  for  convenience  of  presentation  only. 


determined  by  a  single  measurement,  the  difference  between 
the  Ne+  and  Ar+  rate  coefficients  for  each  anion,  while  all 
other  reported  values  benefit  from  redundant  measurements. 

DISCUSSION 

This  experimental  study  has  produced  a  set  of  12  rate 
coefficients  for  thermal  reactions  of  rare-gas  cations  with 
atomic  halide  anions  at  room  temperature,  comprising  both 
MN  and  TI  channels.  The  reactions  with  Br2“  are  a  byprod¬ 
uct  of  this  work  and  will  not  be  discussed  further,  except  to 
say  that  the  measured  values  with  Ar+  and  Kr+  are  in  agree¬ 
ment  with  the  literature  data.®’  There  are  several  clear  trends 
in  the  data  (see  Fig.  4(a)):  The  reactions  of  the  lighter  rare  gas 
cations  are  faster  than  those  of  the  heavier  rare  gas  ions,  while 
reactions  involving  the  light  halide  anions  tend  to  be  slower. 
The  rate  coefficients  across  the  12  reactions  differ  by  more 
than  a  factor  of  at  least  ~60,  with  some  of  the  reactions  so 


TABLE  III.  Derived  total  rate  coefficients  for  the  indicated  Rg^  -|-  X~  re¬ 
actions  (in  10^*  cm^  s“*  units)  with  the  transfer  ionization  rate  coefficient  in 
parentheses. 


Rg+ 

X 

cr 

Br- 

r 

Br2 

Ne+ 

3.0  ±  l.D 

3.3  ±  1.2““ 

6.5  ±  1.7“ 

Ar+ 

0.4  ±  0.2’’ 

0.5  ±0.2  (0.19) 

2.1  ±0.4  (0.11) 

3.7  ±  0.8 

Kr+ 

<0.2’’ 

<0.2’’ 

0.9  ±  0.2“ 

3.2  ±0.8 

Xe+ 

<0.1’’’" 

<0.3’’ 

<0.2’’ 

4.0  ±  1.5 

^Transfer  ionization  is  exothermic,  but  the  rate  coefficients  have  not  been  determined. 
’^Transfer  ionization  is  endothermic. 

‘^The  rate  coefficient  for  Xe”*"  +  Cl“  was  fixed  in  the  analysis  to  a  value  <10“^  cm^  s“^ ; 
see  text. 
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slow  that  only  upper  limits  could  be  determined.  In  contrast, 
reported  rate  coefficients  for  polyatomic-monatomic  MN 
reactions  differ  by  no  more  than  about  a  factor  of  3.® 

These  findings  are  compatible  with  accidentally  resonant 
charge  transfer  in  which  the  attached  electron  is  transferred  to 
excited,  nearly  resonant  states  of  the  rare  gas  atom.  Suitable 
excited  states  of  the  rare-gas  atoms  are  their  metastable  states 
and  those  above  because  they  lie  below  the  Rg  ionization  limit 
by  approximately  the  electron  affinity  of  the  halogen  atoms 
(~3-3.6  eV).  Similar  arguments  are  invoked  in  the  context 
of  thermal-energy  reactions  of  doubly-charged  rare-gas  ions 
with  rare-gas  atoms. This  suggests  that  the  Rg+  -|-  X“ 
reaction  rates  should  be  sensitive  to  the  energy  difference 

AE  =  lE(Rg)  -  E*(Rg)  -  EA(X),  (4) 

where  lE(Rg)  is  the  ionization  energy  of  Rg,  E*(Rg)  is  the 
energy  of  the  lowest  excited  state  of  Rg,  and  EA{X)  is  the 
electron  affinity  of  X.  According  to  Olson,^®  Landau-Zener 
transitions  in  MN  systems  are  most  effective  when  the  initial 
and  final  states  intersect  at  crossing  distances  R^  from  15  to 
50  ao.  Because  the  initial  Rg+  -|-  X“  curve  is  dominated  by 
the  Coulomb  potential  and  is  essentially  system  independent 
at  large  distances,  Rj^  may  be  calculated  from  only  the  asymp¬ 
totic  energy  difference  between  the  initial  and  final  states,  AE. 
In  units  of  Uo  for  R^  and  eV  for  AE,  R^  —  27.2/AE,  mean¬ 
ing  that  Olson’s  favorable  crossing  distances  are  equivalent  to 
AE  —  0.5-1. 8  eV.  Our  experimental  data  (see  Fig.  5)  agree 
with  this  qualitative  expectation:  The  measured  rate  coeffi¬ 
cients  fall  off  rapidly  at  crossing  radii  above  20  ao.  Note  that 
the  estimated  crossing  radii  refer  to  the  lowest  excited  Rg* 
states;  those  for  the  actual  dissociating  states  may  be  larger. 

Hickman'®  derived  a  simple  formula  to  estimate  rate 
coefficients  for  MN  reactions,  and  he  shows  that  it  works  well 
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FIG.  5.  Rg"*'  -I-  X~  rate  coefficients  from  experiment  (closed  circles)  and  es¬ 
timated  using  Hickman’s  parameterization,  Eq.  (5)  (open  squares)  as  a  func¬ 
tion  of  the  crossing  radius  Rj. 


in  many  cases.  His  formula, 

kMN  =  5.3  X  10-^  ^  (5) 

yielding  k  (in  cm^  s“'  units)  has  the  virtue  that  it  depends  only 
on  easily  available  parameters,  the  temperature  T  (in  Kelvin), 
the  reduced  mass  /z  of  the  reactants,  and  the  EA  of  the  halide 
(in  eV).  The  equation  with  slightly  modihed  values  was  found 
to  describe  well  the  variation  in  MN  rate  coefficients  of  poly¬ 
atomic  anions  with  monatomic  cations.®  However,  the  deriva¬ 
tion  assumes  that  the  initial  ionic  potential  curve  crosses 


many  hnal  neutral  curves  and  that  some  of  them  will  occur  at 
favorable  crossing  radii  Rx,  an  assumption  that  breaks  down 
for  some  of  the  reactions  here.  For  instance,  the  Xe+  -|-  Cl“ 
reaction  is  exothermic  by  about  8.5  eV,  but  only  a  handful 
of  product  curves  encompassing  3  electronic  states  of  Xe 
and  2  of  Cl  are  energetically  accessible.  As  may  be  seen  in 
Fig.  5,  the  estimating  formula  comes  within  a  factor  of  ~2 
to  the  measured  values  for  systems  which  happen  to  have  a 
favorable  Rx  ~20  ao,  but  seriously  overestimates  the  rate  co¬ 
efficients  at  larger  Rx,  those  that  correspond  to  small  AE.  If 
one  is  aware  of  this  limitation,  the  formula  continues  to  be 
useful.  Indeed  the  exceptionally  strong  correlation  between 
rate  coefficient  and  Rx  seen  in  Fig.  5  suggests  that  crossing 
radius  is  the  dominant  factor  contributing  to  the  likelihood  of 
MN.  Polyatomic  systems  are  almost  assured  to  have  a  cross¬ 
ing  in  the  favorable  region,  resulting  in  the  observed  small 
variability  between  MN  rate  coefficients. 

A  number  of  MN  rate  coefficients  for  reactions  of 
diatomic  and  polyatomic  anions  with  monatomic  cations  were 
previously  reported  as  a  function  of  temperature  from  300  to 
550  K  using  the  VENDAMS  technique.®  All  of  the  rate  co¬ 
efficients  were  measured  via  the  difference  with  either  Ar+ 
or  Kr+  with  Cl~  or  Br“ .  At  the  time  of  that  study,  those  rate 
coefficients  for  those  reactions  were  not  known,  and  the  data 
were  analyzed  under  the  assumption  that  the  rate  coefficients 
were  less  than  the  uncertainty  of  the  experiment  of  about 
2  X  10“®  cm^  s“'.  While  this  appears  to  be  true  in  the  re¬ 
actions  involving  Kr+,  the  Ar+  reactions  occur  slightly  faster. 
In  light  of  the  current  results,  the  set  of  MN  rate  coefficients 
reported  in  that  study  involving  Ar+  should  be  uniformly 
increased  by  4  x  10“®  cm^  s“^. 

It  has  also  been  possible  to  determine  product  branching 
ratios  between  MN  and  TI  for  some  but  not  all  reactions.  For 
about  half  of  the  systems  studied,  TI  is  endoergic  (Table  111); 
for  several  other  systems  background  X+  prevented  a 


FIG.  6.  Schematic  of  potential  energy  curves  relevant  to  MN  and  TI  relative  to  Ar  +  X  for  the  X  =  I  and  X  =  Br  systems.  Ar*  refers  to  the  lowest  metastable 
state  and  the  electron  energy  in  the  TI  final  state  is  taken  to  be  zero. 
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determination  of  the  TI  rate  coefficient.  A  measurable  contri¬ 
bution  of  TI  has  been  found  for  the  reactions  of  Ar"*"  with  Br” 
and  1“.'^  While  the  two  TI  rate  coefficients  differ  by  less  than 
a  factor  of  two,  the  branching  fractions  varies  much  more; 
TI  accounts  for  ~40%  of  Ar+  -|-  Br“ ,  but  just  ~6%  for  Ar’*’ 
-t- 1“ .  Although  the  dynamics  of  both  TI  and  MN  are  unclear, 
it  is  worth  noting  that  the  smaller  TI  branching  fraction  in  the 
latter  reaction  may  be  a  consequence  of  the  differing  energies 
of  the  final  states. 

Figure  6  shows  approximate  potential  energy  curves  of 
the  relevant  ArX  (X  =  Br,I)  states  relative  to  the  separated  Ar 
and  X  atoms.  The  potential  curve  leading  to  TI  lies  ~1.2  eV 
below  that  for  MN  in  the  I~  reaction,  but  ~0.6  eV  above  the 
MN  curve  in  the  Br-reaction.  Hence,  a  colliding  Ar+  /I“  pair 
will  first  cross  the  MN  curve,  while  an  Ar+/Br”  pair  will  first 
cross  the  TI  curve.  If  one  now  makes  the  assumption  that  the 
autoionization  leading  to  TI  is  most  efficient  when  the  ejected 
electron  has  zero  energy  (plausible  as  the  inverse  process 
e  +  At  +  X+  ^  At+  +  X"  is  a  recombination  reaction 
known  to  have  a  maximum  cross  section  at  low  electron  ener¬ 
gies),  then  TI  will  occur  mainly  at  the  crossings  between  the 
ionic  and  TI  curves  and  it  would  follow  that  TI  is  less  impor¬ 
tant  in  the  I~  reaction  because  the  MN  reaction  reduces  the 
flux  reaching  the  TI  crossing. 

CONCLUSION 

Mutual  neutralization  rate  coefficients  (Table  III)  have 
been  derived  for  reactions  of  rare  gas  cations  (Ne+,  Ar+, 
Kr+,  Xe+^  with  halide  anions  (Cl~,  Br“,  H)  for  the  first  time. 
Unlike  MN  reactions  involving  polyatomic  species,  which 
vary  by  no  more  than  a  factor  of  3,  the  rate  coefficients  for 
the  monatomic-monatomic  processes  vary  by  at  least  a  factor 
of  60.  Rate  coefficients  increase  down  the  periodic  table 
(Cl“  to  I“)  for  the  halide  anions,  and  up  the  periodic  table 
(Xe+  to  Ne+)  for  the  rare  gas  cations.  Previously  reported 
parameterizations  fail  to  predict  the  rate  coefficients  of  these 
reactions,  likely  due  to  the  breakdown  of  the  assumption  that 
the  reactions  involve  a  large  density  of  neutral  product  states. 
Instead,  the  rate  coefficients  vary  strongly  according  to  the 
exothermicity  corrected  by  the  energy  of  the  lowest  excited 
state  of  the  rare  gas  neutral;  systems  with  such  an  adjusted 
exothermicity  of  around  1  eV,  corresponding  to  a  crossing 
of  the  neutral  and  ion  potential  curves  of  ~20  Uq,  show  a 
large  rate  coefficient  comparable  to  those  for  MN  reactions 
involving  polyatomics.  For  several  of  the  reactions,  compe¬ 
tition  between  mutual  neutralization  and  transfer  ionization 
is  observed.  For  those  reactions  where  it  is  possible  to  derive 
the  TI  rate  coefficient,  only  the  reaction  of  Ar+  with  Br“ 
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has  an  appreciable  fraction  of  the  reaction  occurring  by  that 
channel. 
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